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ABSTRACT 

We present high resolution imaging of the low order ( J = 1 and 2) CO line emission from the z = 3.93 
submillimeter galaxy (SMG) MM18423+5938 using the Expanded Very Large Array, and optical 
and near-IR imaging using the Canada-France-Hawaii Telescope. This SMG with a spectroscopic 
redshift was thought to be gravitationally lensed given its enormous apparent brightness. We find 
that the CO emission is consistent with a complete Einstein ring with a major axis diameter of 
~ 1.4", indicative of lensing. We have also identified the lensing galaxy as a very red elliptical 
coincident with the geometric center of the ring and estimated its photometric redshift z ~ 1.1. A 
first estimate of the lens magnification factor is m ~ 12. The luminosity Lco(i-o) of the CO(l — 0) 
emission is 2.71 ± 0.38 x 10 11 mT 1 K km s _1 pc 2 , and, adopting the commonly used conversion 
factor for ULIRGs, the molecular gas mass is M(Ha) = 2.2 x 10 11 mT 1 Mq, comparable to unlensed 
SMGs if corrected by m ~ 12. Our revised estimate of the far-IR luminosity of MM18423+5938 is 



2 x 10 



13 



< L Fm < 3 x 10 1 



L Q , comparable to that of ULIRGs. Further observations are 



required to quantify the star formation rate in MM18423+5938 and to constrain the mass model of 
the lens in more detail. 

Subject headings: submillimeter: galaxies - gravitational lensing: strong - galaxies: evolution - galax- 
ies: starburst - galaxies: ISM 



1. INTRODUCTION 

There is mounting evidence that large ellipical galax- 
ies form the majority of their stars at early epoc hs 
(z > 1), and quickly (duration < 1 Gvn lRenzinil ((2006)). 
Bright submillimeter-selected galaxies (SMGs; defined 
as Ss50fim > 3mJy) are an important galaxy popula- 
tion in this regard, likely representing a key star forma- 
tion epoch with extreme starbursts thought to be driven 
by major mergers of gas-rich galaxies at high redshif t 
(|Smail et al.l2002l:lBlaui^t^[200irTacconi et al1[2U08h . 
T he discovery of evolved galaxies at h i gh redshift (z > 
2: IDaddi et aJ J 120051 iKriek et al.ll200l iToft et all 120091: 
van Dokkum et all 120081: Kurk et all 120091: iCollins et all 
2009t iWilliams et all 120091 iKotilainen et all I2009D . as 
w ell as luminous SMGs at even higher re dshifts (z > 
5: IRiechers et al.l I201CH : ICapak et al.l i2011) provides evi- 
dence for this early active star formation phase out to 
about 1 Gyr after the big bang. 

The very bright SMG MM18423+5938 (S^oo^m = 
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30 ± 3 mJy) was discovered serendipitously at the 
IRAM 30m tele s cope w ith the bolometer array MAMBO 
(|Lestrade et al.l I2009D . Although no optical coun- 
terpart was known for MM18423+5938, its redshift 
z = 3.92960 ± 0.00013 was spectroscopically measured 
through the detection of several strong CO and CI emis- 
sion lines with the wide spectro meter EMIR also at this 
telescope (jLestrade et all 120100 . These line detections 
indicate a massive reservoir of molecular gas in a galaxy 
only 1.6 Gyr after the big bang. The observed CO ladder 
(CO(4-3), CO(6-5) and CO (7-6), CO(9-8) upper limit) 
peaks around J = 5, indicative of a starburst-powered 
SMG that is not dominantly excited by an active nucleus 
(AGN). 

Recently, bright, strongly lensed SMGs have been 
discovered in large area surveys jVieira et al.l l2010t 
iNegrello et alll2010l: iFraver et al.ll20lH) . Even compared 
to these, MM18423+5938 is very bright, and its well 
established redshift makes it a prime target for high- 
resolution observations to study the cause of its high 
brightness, to understand its place in this emerging 
population, to verify the hypothesis that it is gravi- 
tationally lensed, and to study further the similarities 
of MM18423+5938 with the other very bright SMG 
SMMJ2135-0102 at z = 2.3 (S H5Qum = 106 mJy) which 
has b een shown to be strongly lensed (Swinban k et all 

MM- 

In this Letter we report observations of the CO(l- 
0) and CO(2-l) line emissions in MM18423+5938 to 
0.6" res olution using the E xpanded Very Large Array 
(EVLA; IPerlev et all (|2011D 1 as well as multicolor op- 
tical and near-IR imaging using the Canada-France- 
Hawaii Telescope (CFHT). These observations confirm 
that MM18423+5938 is lensed, providing a unique op- 
portunity to study the gas dynamics on small physical 
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Fig. 1. — EVLA observations. Top left : spectrum of the CO(l-O) emission from MM18423+5938 at 51 km s -1 resolution. Bottom left : 
spectrum of the CO(2-l) emission at 26 km s -1 resolution. Zero velocity corresponds to the redshift z = 3.9296. Right : image of the 
CO(2-l) emission from MM18423+5938 integrated over 230 km s -1 at 0.6" resolution. The rms is 0.16 mjy beam -1 , and the contour 
levels are -0.4,-0.2, 0.2, 0.4, 0.6, 0.8, 1.0, 1.2 mjy beam -1 . Negative contours are dashed. The cross shows the position of the lensing 
galaxy at the ring center. 



TABLE 1 

Observed line parameters of 2-component Gaussian fit 



line S„ AVfwhm Vo Ico l 'co 10 " 

(mjy) (km/s) (km/s) (Jy km/s) (K km/s pc 2 ) 

CO(1-0) 2.05 ±0.22 161 ±30 -14 ±11 0.35 ± 0.05 2.71 ± 0.38 
0.62 ±0.38 118 ± 114 -228 ± 52 0.078 ± 0.06 

CO(2-l) 10.83 ±0.74 239 ± 26 -6 ± 9 2.74 ±0.23 4.74 ± 0.40 
1.92 ± 1.34 125 ± 128 -286 ± 59 0.25 ± 0.20 

Note. — Zero velocity corresponds to a redshift z = 3.9296. The apparent 
L'qq results from the two components added up. Only statistical uncertainties are 
given. 



scale with further high sensitivity and high spatial res- 
olution EVLA observations. We use a concordance, flat 
A-CDM cosmology throughout, with Ho — 71 km s -1 
Mpc -1 , Q M = 0.27, and A = 0.73. 

2. OBSERVATIONAL DATA AND REDUCTIONS 

2.1. Expanded Very Large Array 

Observations were made in the C configuration of the 
EVLA in December 2010. We targeted the CO(1-0) line 
(rest frequency = 115.271 GHz) at 23.38 GHz, and the 
CO(2-l) line (rest frequency = 230.538 GHz) at 46.76 
GHz. A total bandwidth of 248 MHz was employed, with 
124 spectral channels. Spectral imaging of the broad line 
of MM18423+5938 at high spectral resolution was made 
possible with the ne w capability of th e WIDAR corre- 
lator of the EVLA (|Perlev et all 120111 ). The total ob- 
serving time was 2 hours at 23 GHz and 8 hours at 47 
GHz. Dynamic scheduling was employed to ensure good 
atmospheric observing conditions. The phase calibra- 
tor w as J1849+6705 and the sw itching cycle was 5 min- 
utes (|Carilfi fc Holdawavl 119991 ) . The synthesized beam 
FWHM using robust weighting (R=l) in the AIPS data 
analysis software was 2.1" x 1.0" at 23 GHz, and roughly 
circular with FWHM = 0.6" at 47 GHz. The flux density 



scale was set by observing 3C48, but some observations 
at 47 GHz were made with the target source and 3C48 
at very different elevations. We estimate an absolute cal- 
ibration uncertainty at 47 GHz of ±20%, even with the 
nominal opac ity correction made using the weather data 
(|Butlerll2010h . The rms noise per channel is 0.16 mjy 
beam -1 at 23 GHz (rest frequency = 115.271 GHz) at 
51 km s -1 resolution, and 0.48 mjy beam -1 at 47 GHz 
(rest frequency = 230.538 GHz) at 26 km s -1 resolution. 

Figure Q] (left panel) shows the CO(1-0) and CO(2-l) 
spectra. Both lines show a possible non-Gaussian tail 
towards lower velocities, so two components were fitted 
to both the CO(1-0) and CO(2-l) profiles. Line peak 
flux densities, line widths, mean velocities and velocity- 
integrated line fluxes Ico are in Table [1] 

Figure Q] (right panel) shows the CO(2-l) emission in- 
tegrated over the channels of the full line width (230 
km s -1 ) at 0.6" resolution. The morphology is consis- 
tent with an Einstein ring interpretation. The ring is 
elliptical, with a major axis ~ 1.4" and a minor axis 
~ 1.2". The J2000 coordinates of the ring center (geo- 
metric center of the inner hole) are a= 18h 42m 22.29s, 
S = 59° 38' 29.5" ± 0.2" in the International Celes- 
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Fig. 2. — CFHT images. RGB color images of 60" X 60" regions centered on MM18423+5938 in the stacked MegaCam g,r, and i, and 
WIRCam J, H , and Kg band images, respectively, are shown in the left and right panels. The pixel scales are 0".187 and 0".304 in the 
MegaCam and WIRCam images, respectively. The main lens galaxy is circled in red and is at the center of the insets (8" X 8"). The 
red dots in the insets are the FWHM point-spread function (PSF) during observations (in r-band, and Ks band). The optical image was 
affected by scattered light from two nearby ~ 9 Vmag stars, see §2.2 for details. 



TABLE 2 

CFHT GRIZ AND JHKg BAND PHOTOMETRY (EXTINCTION CORRECTED AB-MAGNITUDES) 



Object g r i 

Lens Galaxy 26.33 ± 0.43 24.84 ±0.20 

Nearby Gal. (N.E.) 23.75 ± 0.05 22.28 ±0.03 

tial Reference System (ICRS). No continuum emission is 
detected at either frequency, with 3a limits of 75 /iJy 
beam" 1 at 23 GHz and 210 /zJy beam" 1 at 47 GHz. 

2.2. CFHT Optical and Near-IR imaging 

To identify the putative lensing galaxy, we carried out 
multi-band optical and Near-IR (NIR) imaging at the 
CFHT with the prime focus imagers, MegaCam and 
WIRCam respectively. The observations were com- 
pleted in 2010 November, with total on-source integra- 
tion times of 4.75h on MegaCam (6300s in i, 5400s in 
r, and 2700s each in g and z filters), and 0.58h on 
WIRCam (950s in K s , 840s in J, and 285s in bands). 
The median seeing during the optical imaging was 0" .65 
in the r-band, which was well sampled by the Mega- 
Cam pixel scale of 0".187; during the NIR imaging, the 
median seeing was 0".7 in Kg-band, marginally sam- 
pled by the WIRCam 0".304 pixel scale. Standard 
calibrations were obtained and the data processed at 
CFHT with the Elixir and I'iwiv.2 pipelines. The sky- 
subtracted dit hered exposures in each filter were s tacked 
using Swarp (jBertin et al.l 12002) and MegaPipe (jGwyrJ 
2008). For obtaining consistent photometry in all pass- 
bands, the NIR images were re-convolved from 0".304 
pixel scale to the MegaCam resolution of 0".187 dur- 
ing stacking of the dithered WIRCam images. Astro- 
metric and photometric calibratio ns in the op tical and 
NIR were c arried out with Sc amp (|Bertinll2006l ). against 
USNO-B1 (IMonet et al.l [2003h. and 2MASS Point Source 
(jSkrutskie et all 12006( 1 catalogs. The 2MASS positions 
are tied to the ICRS system via the Hipparcos Tycho 
reference frame to better than 0.2" (2MASS Explane- 
tary Supplement Chap I, subsection 6b, paragraph ix). 
There are thirteen 2MASS sources in our full WIRCam 



z J H K s 

21.81 ±0.24 21.25 ±0.33 20.74 ± 0.08 
20.67 ±0.04 19.19 ±0.03 19.14 ±0.04 19.08 ± 0.02 



image that have been used for this NIR-radio registra- 
tion. The RGB color stamp images (Figures [2]) con- 
structed with optical gri filters, and the NIR JHKs 
bands, show a red galaxy (with ellipticity = 0.19) at 
a=18h 42m 22.27s±0.29", 5 = 59° 38' 29.6" ± 0.25" in 
the ICRS, i.e. in the same coordinate system as our 
EVLA map. In the NIR (Fig right), the detection sig- 
nificance is at 10er in Kg, and at ~ 3c in J and H bands. 
In the optical (Fig[2J left), the bright stellar halo from 
the two nearby ~ 9 Vmag saturated stars led to signifi- 
cant background noise, and therefore to lower detection 
significances but nonetheless successful identifications in 
the i,g,r filters. The z filter was affected by fringing. 

Given the low detection significances in the optical 
imaging, magnitudes we re measured using S extractor 
(jBertin fe Arnoutslll996f) in dual image mode, using the 
Kg band image for object detection, while measuring the 
flux in matched apertures in all filters; for the optical fil- 
ters, the WIRCam Kg image was resampled to match the 
pixel scale of MegaCam. Further, we carried out photom- 
etry of the nearby galaxy at a — 18h 42m 22.52s± 0.36", 
S = 59° 38' 33.17" ± 0.31" in the ICRS (lying 5" E of 
N of the putative lens, see Fig [2]) expected to contribute 
to the lensing model. The foreground Galactic extinc- 
tion corrected ^45-magnitudes are in Table [2] (extinction 
E(B - V) is 0.0462 ± 0.0009 at the source position in 
Schleger maps). 

3. ANALYSIS AND DISCUSSIONS 

Figure [3] shows an overlay of the CO(2-l) contour map 
of MM18423±5938 and the Kg band image. The main 
lens galaxy is coincident with the geometric center of the 
ring. With the measured NIR a nd optical magnitud es, 
we used SED fitting in LePhare (|Arnouts et al.lll999| ) to 



24.17 ±0.15 
21.33 ±0.02 
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estimate a photometric redshift z p hot ~ 1-1, with a for- 
mal uncertainty of ±0.1, for the lens best fit with an 
early type galaxy template. No satisfactory fit was ob- 
tained for the North-East nearby galaxy mainly because 
the J magnitude did not match with the other bands. 
We built a preliminary lens model based on isother- 
mal elliptical potential s for the two galax ie s using the 
LENSTOOL Software (|Kneib et al.1 lfl996ft . Uullo et al.l 
(|2007t) V This resulted in a quadruple image with the 
two brightest spots located at the two maxima of the 
CO (2-1) map, and yielded the total magnification factor 
m ~ 12. The nearby galaxy contributes to the lensing 
model in making the ring slightly ellipitical as observed 
but its mass and redshift cannot be separated with our 
current data. Had we derived its redshift, we could have 
speculated whether or not the J band excess could be 
due to strong emission of its H a line. The physical scale 
probed in the source with the beam of our EVLA map is 
0.6" x (7.1kpc arcsec -1 )/-^/^, i.e. 1.2 kpc with our first 
estimate of the magnification factor m. 

We revise the far- IR luminosity Lwi r of 
MM18423±5938 given in ILestrade etaLl (|20iri since 
the strong Spitzer MIPS 70/im source, tentatively 
associated with MM18423±5938 then, has now a 
position that is discrepant by more than 5er with 
the more accurate EVLA position newly found for 
MM18423+5938. The SED of the dust emission was 
modeled by a single temperature modified black-body 
function with the standard power law for the dust 
opacity, fitted to the av ailable photome t ry at A= 1.2, 
2 and 3 mm given by ILestrade et al.l (|2010f ) and a 
3-cr upper limit of 3 mJy at 70/im. We varied the 
opacity index over the range 1.0 < j3 < 2.0 and the 
temperature ov er 20 < Tnnst < 60 K that are plausible 
for SMGs (e . g. iMagnelli et al.l 12015 llvisqn et al.l l20Tot 
iHwang et all 120101: iChapman et al.l [2Q10), and found 
the range 2 x 10 13 mT 1 < L FJR < 3 x 10 14 to -1 L 
with a 95% confidence level, where Lfir is mostly 
the luminosity between A ~ 70/im and 3mm, and m 
is the magnification factor. Using the FIR to star 
formation rate conversion from iKennicuttl Jl998), this 
luminosity corresponds to a range in star formation 
rate of 3 300 m" 1 < SFR < 22000 m" 1 M^yr" 1 . 
Recently, M cKean et al.l (|2011| ) made a deep image 
at 1.4 GHz with the Westerbork array (WSRT) and 
detected MM18423±5938 at a position consistent with 
ours to within la. They used the radio-FIR relation to 
infer an apparent luminosity Lfir — 5.6^24 x 10 13 L© 
consistent with our new estimate. 

The velocity-integrated CO(1-0) flux density yields 
a luminosity L' co(1 _ 0) = 2.71 ± 0.38 x 10 11 m" 1 K 

km_s_~ 1 pc 2 using the relationships from Solomon et al. 
(|1992f) . If corrected by m ~ 12 this is the typica l 
luminosity f o und f or SMGs, e.g. , iHarris et al.l (|2010[) . 
iFraver et all (|201lD . Ilvison et all ((20111). This implies 
an H2 mass of 2.2 x 10 11 m" 1 M© in assuming a 
CO(1-0) luminosity to H2 mass conversion factor of 
a = 0.8 M Q / K km s" 1 pc 2 , s uitable for starburst 
galax ies in the nearby Universe (iDownes &: Solomon 
119981) and typically used for SMGs. e.g. iTacconi et a" 
(2008). The gas-to-dynamical mass ratio exceeds unity 
for the dynamical mass (4.8 x 10 10 mT 1 M©) computed 
for an hypothetical edge-on disk by using the FWHM of 




26"- 



v . ■ ■ ■ T.-.i 

18 h 42 m 22.75 s 22.50 s 22.25 s 22.00 s 21.75 s 

Right Ascension (J2000) 

Fig. 3.— Overlay of the CO(2-l) contour map of MM18423+5938 
and the Ks band image of the main lens galaxy and N.E. nearby 
galaxy (gray scale). Contours are the same as in Fig[l](right panel). 
Pixel size of the Ks band image is 0.304". 

the lines (~ 200 km/s) as an indication of the rotation 
velocity. This exceedingly large ratio and the fact that a 
disk in MM18423+5938 is not expec ted to be supported 
by the radiation pressure of an AGN (|Scoville et al.|[l995l) 
suggest that a disk would be inclined. Note that this 
ratio still exceeds unity if a merger model is used in- 
stead, the enclosed dynamical mass is approximat ely a 
factor of 2 larger in this case (|Genzel et al.ll200l . In 
Fig[TJ the second component in the CO(1-0) and CO(2- 
1) lines detected at a low significance level is blueshifted 
by 200-300km s _1 . We have verified that this com- 
ponent is not visible in the averaged profile of the three 
high order (J=4 6,7) C O lines observed at the IRAM 30m 
(|Lestrade et al.|[2010D . which indicates there are distinct 
CO components that have low excitation only within 
MM18423+5938. Such a component would be indica- 
tive of either a more extended low order CO component 
of the rotating disk or possibly the second component of 
a merger. 

The high order CO excitation ladder observed with 
EMIR/IRAM was fitted with a large velocity gradient 
(LVG radiative transfer model, T/. = 45K, n(H2) = 
10 3 cm- 3 , N C o/AV = 3 x 10 18 cm- 2 /(km s" 1 )) in 
Lestrade et al. (2010). The CO(1-0) peak flux den- 
sity (2.1 ± 0.2 mJy) is consistent with this model, while 
the CO(2-l) measurement (10.8 ±2.2 mJy) is marginally 
in excess by 54 ± 30% when the 20% calibration uncer- 
tainty at 47 GHz is included. Although this uncertainty 
is presently large, we note that the ratio of the flux den- 
sities CO(2-1)/CO(1-0) is difficult to explain even with 
differential lensing if the excitation is spatially variable in 
the galaxy, unless CO (1-0) is moderately optically thick, 
with rising optical depth toward CO(2-l). New EVLA 
observations are required to investigate this in more de- 
tail. 

The peak rest-frame brightness temperature of the ring 
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is 12 K (East clump in our map of Fig [1}. This is similar 
to wha t has been seen in other z ~ 4 SMGs, such as 
GN20 (jCarilli et all 120101) . although we cannot rule-out 
higher Tb clumps that are not resolved by our present 
observations. 



4. CONCLUSIONS 

While there is a long history of CO observa- 
tions from strongly lensed, high redsh i ft galaxies (e.g. 
Alloin et aT]|1997t iBarvainis et alj|1994 iSwinbank et al.l 
20101 ) . MM18423+5938 is only the second complete Ein- 
stein ring seen in CO emission , after the z=4.1 q uasar 
host galaxy PSS J2322+1944 (jCarilli et al.ll200l . We 
identified the lens and derived a magnification factor of 
*~ 12 using a preliminary lensing model. We measured 
a luminosity L' co ^ 1 _ j, and corresponding H2 mass, for 

MM18423+5938 that, when corrected for magnification, 
are similar to the unlensed SMGs. However, the rel- 
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atively large range for Lpjn in our present analysis re- 
quires further observations to better sample the FIR SED 
in order to better constrai n the implied s tar fo rmation 
rate. As demonstrated by IRiechers et all ()2008D for the 
molecular Einstein ring observed towards the QSO PSS 
J2322+1944, the brighter SMG MM18423+5938 pro- 
vides a unique opportunity to study the gas dynamics on 
physical scale as small as ~ 100 pc with further high sen- 
sitivity and high spatial resolution EVLA observations. 
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